function y = CL_a(p, x, z)
  
  if x(12) < -p.b  % ground effect
      c = [0 -1236.6 504.81 -72.228 -0.9905 6.3825 0.1453];
      e = [z.^ 6, z.^ 5, z.^ 4, z.^ 3, z.^ 2, z.^ 1, z.^ 0]';
  else
      c = [0 0 0 -65.2 -0.9869 7.084 0.152];
      e = [z.^ 6, z.^ 5, z.^ 4, z.^ 3, z.^ 2, z.^ 1, z.^ 0]';
  end
  
  y = c * e;
  
end